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Abstract 

We investigate the mobility of nanometer-size solutes in water in a uniform external electric field. General arguments 
are presented to show that a closed surface cutting a volume from a polar liquid will carry an effective non-zero surface 
charge density when preferential orientation of dipoles exists in the interface. This effective charge will experience 
a non-vanishing drag in an external electric field even in the absence of free charge carriers. Numerical simulations 
of model solutes are used to estimate the magnitude of the surface charge density. We find it to be comparable 
to the values typically reported from the mobility measurements. Hydrated ions can potentially carry a significant 
excess of the effective charge due to over-polarization of the interface. As a result, the electrokinetic charge can 
significantly deviate from the physical charge of free charge carriers. We propose to test the model by manipulating 
the polarizability of hydrated semiconductor nanoparticles with light. The inversion of the mobility direction can be 
achieved by photoexcitation, which increases the nanoparticle polarizability and leads to an inversion of the dipolar 
orientations of water molecules in the interface. 
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1. Introduction 

Electrophoretic mobility is the drag experienced by 
a dissolved, usually colloidal, particle in a uniform ex¬ 
ternal electric field. Mobility of oil drops and air bub¬ 
bles in water has been known for a long time jT|] and 
is traditionally linked to preferential adsorption of ions. 
Their counterions form the diffuse double layer. The 
overall charge measured by mobility is determined by 
an incomplete compensation between the charge of the 
adsorbed ions and the part of the diffuse layer within 
the shear surface. The latter encircles the stagnant layer 
of the electrolyte moving together with the dissolved 
particle. While the overall force acting on the ions of 
the electrolyte is zero, the electrokinetic drag is the re¬ 
sult of choosing a limited volume within the electrolyte, 
surrounding the colloidal particle, with an uncompen¬ 
sated charge. The dragging force is thus the product of 
the average charge (Qr) within the shear surface with 
the electric field acting on the charges. We show here 
that the idea of a limited volume cut from the liquid 
and producing an excess charge can be extended to the 
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dipolar polarization of the interface. While the dielec¬ 
tric surrounding the nanoparticle is neutral overall, like 
the electrolyte in the standard models, the divergent po¬ 
larization of the interface produces an uncompensated 
bound charge when integrated over a finite volume. 


The excess of the adsorbed charge over the diffuse- 
layer charge, i.e., uncompensated charge (Q)r + 0, is 
reflected in the sign of the ^"-potential at the shear sur¬ 
face §. A negative ^"-potential, typically recorded for 
oil drops and air bubbles in water, has been attributed 
to the excess of the adsorbed negative charge, with the 
hydroxide anion being a long-time favorite yj, |4j, |5j, |6[]. 


Recent calculations 11 i and measurements by 
surface-sensitive second-harmonic generation tech¬ 
niques §111111 m do not support excessive adsorp¬ 
tion of hydroxide to the oil-water 110] and air-water 
10 interfaces. In addition, the total X-ray reflec¬ 
tion fluorescence spectroscopy Cl provides the upper 
estimate for the free surface charges at the air-water in¬ 
terface at the level of 0.002 (e/nm 2 ). Depending on the 
pH and other conditions, this estimate is up to two or¬ 
ders of magnitude below the surface charge density of 
0.02-0.4 (e/nm 2 ) extracted from mobility rail ToL ll . 
It seems plausible that either the formalism of estimat¬ 
ing the surface charge density from mobility requires 
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modification or alternative mechanisms of mobility, not 
involving ion adsorption, might be involved. 

The possibility of charge-free electrophoretic mo¬ 
bility in water has been discussed in the literature 
0 B H d G3 El 0. The main idea here is 
that the microscopic structure of the interface, allow¬ 
ing molecular order within the hydration layers, can ei¬ 
ther produce an effective electrokinetic charge, not re¬ 
lated to charges of free carriers, or substantially mod¬ 
ify the effect of adsorbed ions on the overall mobility. 
This proposal has faced two difficulties. From the the¬ 
oretical side, there is no established framework of how 
to translate the microscopic structure of the interface, 
captured by atomistic numerical simulations, into the 
macroscopic current. Care is required in implement¬ 
ing correct cutoff/boundary conditions |2C)j, 21, 221 and 
statistical ensembles adequately representing the con¬ 
ditions of mobility measurements (as discussed briefly 
below). In addition, the field strengths required to pro¬ 
duce sufficient sampling in simulations are significantly 
higher than experimental fields |23jl and can potentially 
modify the structure of the solution. From the experi¬ 
mental side, it is not clear how to connect the results of 
surface-sensitive spectroscopies, which directly report 
on the polarization structure of the interface [24| 125], 
with measured mobilities. 

Here we address the calculation of the force acting on 
a nanometer-size particle dissolved in water and placed 
in a uniform external field. We do not directly calculate 
the current produced in response to the external field as¬ 
suming that, once the force is known, the mobility can 
be calculated by applying standard equations of hydro¬ 
dynamics ||26l) 1 lh (as shown for the capillary flow in the 
Appendix). Mobility of the hydrated solute is typically 
expressed, through Smoluchovski’s equation (Eq. (JSJi 
below), in terms of either the ^-potential or its effective 
charge. We derive a relation between the effective mo¬ 
bility (electrokinetic) charge and the interfacial struc¬ 
ture of the water dipoles represented by the first-order 
orientational order parameter of the interface. This pa¬ 
rameter is in principle accessible by surface-sensitive 
spectroscopies [24, l25Ll27ll28ll and by equilibrium com¬ 


puter simulations of solutions (29 

This model shows that the effective charge of the so¬ 
lute responding to the uniform external field is not equal 
to the charge of the free carriers. It is therefore possi¬ 
ble that the effective electrokinetic charge reported by 
mobility measurements significantly overestimates the 
number of adsorbed ions. The orientational structure 
of interfacial dipoles is the key in understanding these 
differences. Since the interfacial structure and dipolar 
orientations in the interface can be altered by modifying 


the solute/substrate |29[ 3T, Tcj], one gains the means 
to experimentally test both the effect of the interface 
on the effective electrokinetic charge and the hypoth¬ 
esis of charge-free mobility. In particular, we suggest 
that changing the polarizability of a (semiconductor) 
nanoparticle by exciting electron-hole pairs can invert 
the sign of the mobility. This effect is driven by the rela¬ 
tion between the orientation of dipoles in the hydration 
layer with the nanoparticle polarizability f 3~3l manipu¬ 
lated by light | 


2. Interfacial structure and particle mobility 

2.1. General arguments 

We start by considering a single spherical ion with 
the charge q at its center and with the radius a. It is 
placed in a polar liquid with the bulk static dielectric 
constant e s . We will further consider a spherical liquid 
sample with the macroscopic radius L and place the ion 
at its center to simplify the geometry. The instantaneous 
charge density in the sample is 


p = Pi+Pb, (1) 

where p,- = qd (r) and p*(r) = Yjj qfi ( r ~ r j) is the den¬ 
sity of bound charge at a given instantaneous config¬ 
uration of the liquid with the atomic partial charges qj 
located at the coordinates r r Based on charge conserva¬ 
tion 0 , pb — -V • P is expressed in terms of the polar¬ 
ization density field P. No specific approximation, such 
as the dielectric boundary value problem, is assumed 
here. The instantaneous polarization field is given by 
the microscopic expression [35. 361 

P(r) = Ya ( r “ r f)“5 V 'Z Q ^( r “ r j) + ■ ■■■ (2) 

j j 

Here, m ; denotes the molecular dipole, Q ( is the molec¬ 
ular quadrupole (defined according to Ref. |f35tl ), and the 
dots refer to the higher-order multipolar terms. When 
the statistical average over the configurations of the liq¬ 
uid is performed, one arrives at statistically averaged 
scalar and vector fields, (pb ) and (P). 

From Eq. ©I, one can calculate the overall charge 
within a spherical volume S1 R with the radius R sur¬ 
rounding the ion at its center 

(Qr) = f \Pi - V ■ (P)] dr. (3) 

Jn R 

By using the Gauss theorem, integration in Eq. © 
yields 

(Qr) =q+ [S n P„ - S R P R ]. (4) 
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Here, P a and Pr are the average radial projections of the 
polarization density, P r — r ■ (P), r = r/r taken at r = a 
and r = R, respectively. Further, S a = 4;rcr and Sr - 
AnR 2 are the surface areas. The polarization is locally 
proportional to the field in the continuum electrostatics 
and P r = (q/Anr 2 ){\ - e _1 ). Therefore, in this case, 
r 2 P r is independent of r and the two summands in the 
brackets in Eq. (Q} cancel out. One gets (Qr) = q. 

The dielectric sample is overall neutral and one can 
additionally require 


I 


( Pb)dr = 0 , 


(5) 


where the integral is taken over the liquid volume Q be¬ 
tween the spheres r — a and r — L. This relation im¬ 
poses the boundary condition 


a 2 P a = L 2 P l , 


( 6 ) 


which is satisfied for continuum electrostatics. 

The normal, outward to the dielectric component of 
the polarization field cr q = P„ = —P a plays the role 
of the surface charge density of a discontinuous dielec¬ 
tric interface 0411 . While this charge originates from 
a divergent polarization of bound molecular charges, it 
is experimentally observable. To show that, one can 
consider the electrostatic potential created by free and 
bound charges inside or outside of the macroscopic 
sample. It is given as a sum of the electrostatic poten¬ 
tials arising from the free and bound charges [37 ] 


= q - + 


s„ 


r s 


■dS =I-(l-lU 

r \ ej r 


(7) 


where the surface integral is over the surface of the ion 
S a ■ The overall potential <p = q/(re s ) is said to be di¬ 
electrically screened. This physically implies that any 
probe charge placed outside of the dielectric sample will 
sense the combined charge q e s - q/e s , resulting from 
adding the ion charge with the opposite bound charges 
non-uniformly distributed around the ion and producing 
a non-zero divergence V • P. 

We now move to the next step to point out that the 
polarization field of liquid interfaces often shows a be¬ 
havior more c omp lex than P r cc r 2 of continuum elec¬ 
trostatics [38;, 39, 40]. The function P r often displays 
overscreening, which means that it can be much larger 
in the magnitude at the contact with the ion than pre¬ 
dicted by dielectric models. It also shows oscillations 
caused by molecular granularity as it decays to the r~ 2 
asymptote at r —> oo. While the overall neutrality con¬ 
dition © still must hold, the charge (Q)r obtained by 
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Figure 1: Cartoon depicting the radial projection of the microscopic 
polarization density P r (solid line) and its dielectric form oc r ~ 2 
(dashed line). The volume integral of dP,-/dr between surfaces r = a 
and r = R in Eq. (3j can be non-zero, while it always vanishes in the 
dielectric limit. 


integrating in Eq. ([3} over a small volume £1r can be 
nonzero for a function P r (r) = p(r)/r 2 with a generally 
oscillatory p(r) such that p(°°) = 1 (Fig. |T]). 

This simple observation is the basis of our proposed 
alteration of the standard model of ionic mobility un¬ 
der the drag of a uniform electric field. We suggest that 
( Q)r 4- q if the liquid within the shear surface, dragged 
along with the ion, carries some molecular interfacial 
structure affecting the radial distribution of the polariza¬ 
tion density. The effective charge associated with mo¬ 
bility is affected by the distribution of the bound charge 
within the shear surface, in addition to the total charge 
of free carriers. 


2.2. Ionic mobility 

The hydrodynamic mobility of an ion is determined 
by the shear surface of the radius R , which is coarse¬ 
grained to smooth out the details of molecular granular¬ 
ity by averaging out the molecular motions on the time 
short compared to the time-scale of hydrodynamic flow 
(Fig. [2}. Electrostatics suggests that the force acting on 
the ion and its stagnant layer is the product of the aver¬ 
age charge (Qr) within the shear surface with the field 
acting on these charges. This field is the cavity field 
E c 113611 combining the field from external charges with 
the field of the polarized dielectric outside of the shear 
surface 

(F) = (Q r )E c . (8) 


In dielectric theories, the cavity field inside a sphere is 
related to the macroscopic Maxwell field E by the equa¬ 
tion lf36ll 

E c = —(9) 

2e s + e„ 


where e p is the dielectric constant of the particle. 
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Figure 2: Ion with the charge q and the radius a immersed in a po¬ 
lar liquid in the uniform macroscopic (Maxwell) field E. R indicates 
the radius of the shear sphere incorporating the stagnant layer of the 
liquid dragged by the field along with the solute. E c is the cavity 
field of the uniformly polarized liquid created inside the shear sphere. 
The arrows on the opposite sides of the spherical cavity indicate wa¬ 
ter dipoles oriented favorably (left) and unfavorably (right) along the 
external field. The difference in the chemical potential between right 
and left is positive. It creates an osmotic pressure pushing the particle 
in the direction opposite to the field and corresponding to an effective 
negative charge. 


The steady flow of dissolved particles with the speed 
u is reached when the electrostatic drag is counterbal¬ 
anced by hydrodynamic friction, (F) = 6nquR, where // 
is the bulk viscosity. The resulting mobility p = u/E ^ 
(Qr)/( 4nr]R) [e, » e p , e s » 1 in Eq. CD] gives di¬ 
rect access to the total charge (Qr). Smoluchovski’s 
equation, typically used in the literature, re-writes this 
relation in terms of the ^-potential defined as the elec¬ 
trostatic potential at the shear surface ( — {Qr)[((sR) 
I 4lh . The result is the equation for the mobility |2[] 


P = 


EL 

4nq 


00 ) 


This formalism is well established, and the results of 
mobility measurements are often cast in terms of the ef¬ 
fective surface charge density cr eff = (Qr)/S , where S 
is the surface area of the particle. We follow this estab¬ 
lished practice and focus mostly on (Qr) and the cor¬ 
responding cr eff . The arguments given here need to be 
modified with the account for the diffuse potential when 
electrolyte is present [41]. We do not consider these ef¬ 
fects here and focus instead on charged or uncharged 
particles dissolved in a polar molecular solvent, which 
establishes a microscopic multipolar structure in the in¬ 
terface. The main outcome of this perspective is the 
modification of the effective charge (Qr) by the dipolar 
order of the interface expressed in terms of the average 
cosine of the interfacial dipoles (order parameter) p\. 

Starting from Eq. Q, one can proceed in two steps 
and first apply a reasonable approximation to the surface 
charge density at the shear surface. The surface charge 
density at the actual physical surface of the solute then 
becomes our main focus. Since the shear surface does 
not involve any physical disruption of the liquid struc¬ 
ture, it is reasonable to assume that Pr can be related 


to the field of the ion charge by the rules of continuum 
electrostatics [36] SrPr - q(e s - 1 )/e s . We stress that 
this assumption does not affect the main points of our 
reasoning, as will be clear from the discussion below. 
With the continuum polarization at the shear surface one 
gets in Eq. (01 


(Qr) = — cr q S a . 


( 11 ) 


Since the microscopic susceptibility of the nanome¬ 
ter interface can significantly deviate from_the rules 
of macroscopic continuum electrostatics [38| [42| So], 
cr q = P„ - —P a is left unspecified in Eq. ( [111 . The 
simple message delivered by Eqs. ( ITOb and ( fTTI) is that 
asymmetry in the water susceptibility between the shear 
and solute dividing surfaces leads to a modification of 
the standard result (Qr) = q. 

Since cr q is given by the normal projection of the 
polarization density in the interface, Eq. ( flTb offers a 
new result typically absent in standard dielectric mod¬ 
els. Those suggest that cr q is proportional to the ion 
charge q. However, if the polar liquid is spontaneously 
polarized in the interface, i.e., if the interfacial dipoles 
possess preferential non-random orientations caused by 
the interfacial order i 39i 43,44], (Qr) + 0 even at g = 0. 
What is required is a nonzero radial projection of the 
dipolar polarization density at the solute surface. 

The dipole ordering in the interface can be described 
by the first-order orientational order parameter p\ = 
(rh • f)„, which is the average cosine of the water dipole 
moment projected on the radial direction and calculated 
in a narrow layer at the solute surface r — a 142. 30]. 
The surface charge density can be written in terms of 
the water dipole moment m and the order parameter, 
-cr q - ( mpi/S)(dN s /dr)\ r a . Here, N s - N s (r) is the 
number of water molecules within the shell of the radius 
r > a. By using the definition of the number of water 
molecules in the shell in terms of the solute-solvent ra¬ 
dial distribution function (RDF) go s (r), one can re-write 


■ pmp\G, 


( 12 ) 


where G = gos(a) is the contact value of the solute- 
solvent RDF and p is the number density of bulk wa¬ 
ter. Equation (IT2l) is written for an arbitrary value of q, 
which means that p\G should be calculated in the pres¬ 
ence of the ion charge q\ ctq corresponds to q — 0. 

The value of cr q for large particles can be estimated 
from the a —> oo asymptote for the hard-sphere (HS) so¬ 
lute 145l] Ghs —> PP/p, which results in -cr q — > / 3mp\P , 
where P is the hydrostatic pressure. This gives for water 
at ambient conditions -cr q — 10 _3 /?i(G/GHs)Cf7at m ) 
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e/nm 2 , where G / G hs — exp[-/iA/r w ] defines the affin¬ 
ity of water toward the solute A/r„. beyond simple HS 
packing preferences. 

Equation (ITU establishes the effective charge of a 
closed spherical interface within a polar liquid. Its sign 
is fully defined by the orientational order parameter p \: 
it is negative when the water dipoles preferentially ori¬ 
ent toward the solute/cavity and is positive when they 
point toward the bulk. This equation shows that any 
closed dividing surface, cutting a volume from a polar 
liquid, will be dragged by an external electric field if a 
preferential orientation of dipoles in the interface is es¬ 
tablished. This result is independent of the presence of 
the electrolyte since bound charges are not screened by 
the ions. 

The proposed formalism equally applies to the prob¬ 
lem of a water drop in a nonpolar solvent (oil) I18I1 . 
Equation (fl2l > still defines the surface charge density 
with the convention that the orientational order parame¬ 
ter is calculated by projecting the surface water dipoles 
on the radial direction pointing toward water (inward in 
the case of a drop). To make our assignment clear, the 
water-oil interface with water’s hydrogen pointing to¬ 
ward the oil phase 010 will, according to Eqs. 
(DU and dm produce a negative charge of the water 
drop. 

It is important to note that the electrostatic force lin¬ 
ear in the external field, ( F ) oc E c , assumes an un¬ 
perturbed orientational structure of the interface pro¬ 
jected on the order parameter p\ 04811 . The relaxation 
of the interfacial order in response to the external field 
would represent the interfacial polarizability, which 
contributes to the overall force as a term quadratic in 
the external field. We do not consider the interface po¬ 
larizability here assuming that macroscopic fields used 
in experiment are weak compared to microscopic inter¬ 
facial fields ll39ll and do not significantly alter dipolar 
orientations in the interface. 

The electrolyte is overall neutral and the overall force 
acting on the electrolyte ions is zero F = q>E = 0. 
However, producing current requires work of the exter¬ 
nal source. The power P, or the rate of doing work, is 
related to the current density j 


P = 


f* 


Edr = P e i + N 0 (Qr)uE, 


(13) 


where P e 1 = ( J+ - J-)E is the power required to move 
the electrolyte ions with the overall current of cations 
and anions given as J ± ; No is the number of colloidal 
particles (see Supplemental Material for detail). Equa¬ 
tions (0 and ( 1 1 31 in principle allow a non-zero cur¬ 
rent and power production at q = 0, i.e., for over- 
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Figure 3: Surface charge density of hard-sphere (HS) and Kihara (KH) 
solutes of varying solute size Rq s in SPC/E water B3l (solid points) 
and TIP3P water (open points). The LJ energy eo* in Eq. 013 was 
varied in the simulations: 0.65, 3.7, and 8.0 kJ/mol. The dotted lines 
connect the points to guide the eye. 

all neutral solutes surrounded by a polarized interface. 
This possibility was viewed in Ref. |49| as contradict¬ 
ing to Saxen relations between the streaming potential 
and electro-osmotic current, which are specific forms of 
the Onsager reciprocal relations |26j |. We show in the 
Appendix that the Onsager relations are obeyed in our 
model by the simple fact of being based on the Coulomb 
law applied to both free and bound charges. 

The drag experienced by a closed surface can be 
viewed as a specific form of osmosis 150]. The gradi¬ 
ent of the chemical potential of interfacial waters at the 
opposite sides of the surface is created by the external 
field. It is the consequence of the favorable orientation 
with the field of the molecules on one side of the surface 
compared to the unfavorable orientation on the opposite 
side 15111 (surface arrows in Fig. [2} The chemical po¬ 
tential gradient will result in the osmotic pressure dif¬ 
ference on the opposite sides of the surface as long as 
spontaneous order in the interface persists. This phys¬ 
ical interpretation of non-zero mobility implies that di¬ 
rect numerical simulations of this effect will require the 
pWT ensemble )45, 5j]], keeping the chemical potential 
of water constant. Since these results are presently not 
available, we use more conventional NVT and NPT sim¬ 
ulations of nonpolar and ionic solutes in water to esti¬ 
mate the interfacial charge density cr q in Eq. (fT2l) from 
the computed p\G parameter. 

3. Computer simulations 

We have considered HS and Kihara (KH) solutes dis¬ 
solved in force-field water. The Kihara potential is the 
HS core modified with the Lennard-Jones (LJ) layer at 
its surface [52]. Specifically, the solute-solvent poten¬ 
tial is given as 


(f>o s (r) = 4eoj 


o~o s 


J?Hs) 


, 12 


0~(b \ 6 

f^HS / 


(14) 
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where eos in the LJ energy and <xo v is the distance be¬ 
tween the solute HS core with the radius Rus and wa¬ 
ter’s oxygen. <xo v = 3 A was kept constant in the simu¬ 
lations, while Rus and €q s were varied. 

The molecular dynamics (MD) and Monte Carlo 
(MC) simulations presented here address the question 
of whether the product p\G characterizing water inter¬ 
facing these solutes will produce cr q comparable to ex¬ 
periment. The experimental cr q reported in the litera¬ 
ture are derived from mobility through Smoluchovski’s 
equation [Eq. (0], The details of the simulation 

protocols have been discussed elsewhere |42, 32[ 53] 
and are given in the Supplemental Material. Here we 
focus only on the results. 

Figure [2 shows cr o/e (q — 0, e is the elementary 
charge) from the simulation data changing with the size 
of the HS and KH solutes in TIP3P and SPC/E water 
models I154ll . The size of the KH solute is measured as 
Ros = Rus + o- 0s [Eq. < 1 1 Oi l. It approximates well the 
position of the first peak of the solute-solvent RDF. The 
size of the HS solute Rq s is defined as the distance of 
the closest approach of the water oxygen to the solute. 
It gives the exact position of the RDF’s first peak. 

The sign of the surface charge density is negative for 
both HS and KH solutes, reflecting the preferential ori¬ 
entation of the surface water dipoles into the bulk. In¬ 
creasing the solute-solvent LJ attraction makes the hy¬ 
dration shell denser, as reflected by a higher -cr q. The 
fast drop of -ctq for the HS solute is caused by its partial 
dewetting [55] when Rq s > 5 A. 

The magnitude of cr 0 is somewhat higher than the 
values typically reported from mobility measurements 
(~ -0.04 (e/nm 2 ) for hexadecane in 0.2 mM NaCl at 
pH = 7 |5[]). We estimated the ("-potential for the 
€q s = 0.65 kJ/mol Kihara solute ||42!| (Fig. 0 It 
turned out that (TqRq.. is an approximately linear func¬ 
tion of Ro s (Fig. S3 in Supplemental Material) such that 
£ - 0.026(e/nm)(A’o s //' 1 ) for large Kihara solutes. Ne¬ 
glecting the difference between Rq s and R in this limit 
results in £ ^ 38 mV. This number is not very differ¬ 
ent in magnitude from those typically reported experi¬ 
mentally. For instance, £ ^ -81 + 14 mV was reported 
for xylene droplets in 10“ 5 M NaCl electrolyte at pH = 6 
ft For water at room temperature, the Debye-Hiickel 
length is k - 1 - 3/c 1 ^ 2 A for a single-charge electrolyte 
with the molar concentration c [56]. For the cited exper¬ 
iment, one gets k _ 1 =* 10 3 A and the amount of counte¬ 
rion charge within the stagnant layer of < 1 nm in thick¬ 
ness 115711 can be neglected. The measured ((-potential 
thus reflects the effective electrokinetic surface charge. 
We stress that our solutes are significantly smaller in 
size than oil drops used in the mobility measurements 
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Figure 4: Surface charge density of HS cations (C) and anions (A) 
in TIP3P water. The calculations are done according to Eq. mi : the 
dotted lines connect the points. cr q is negative for cations and positive 
for anions. 


(~ 100 112] to ~ 200 - 300 liToh nm) and have a smooth 
surface, in contrast to the corrugated surface of oil emul¬ 
sions |47|. 

The experimental ("-potential (ll] has the sign oppo¬ 
site to that calculated for the Kihara solutes. The reason 
is the positive sign of p\ in the Kihara-water interface, 
while negative p\ values have been recently reported for 
the oil-water interface Eil. The access to water ori¬ 
entation in the interface is experimentally provided by 
heterodyne-detected vibrational sum-frequency genera¬ 
tion (VSFG) spectroscopy through the imaginary part 
of the VSFG signal Imy (2) ifz/l ujl l. The combination 
of the sign of Im^ (2) and its intensity in principle gives 
access to p\, although in reality fitting of simulations to 
experimental spectra is required IJsIl . Resolving all fea¬ 
tures of the experimentally reported Im^ (2) requires in¬ 
cluding three-body interactions in the force field model 
of water [58]. Whether the same is true regarding the 
values of p\ is not clear at the moment, although there 
are indications that three-site models of water somewhat 
overestimate its spontaneous orientational structure in 
the interface [ 43ll . In addition to spontaneous orienta¬ 
tion in the uncharged interface, the orientation of water 
dipoles and correspondi ng p \ are strongly affected by 
the presence of ions Ii28l,l59i] as we discuss next. 


The potential situation with hydrated ions is illus¬ 
trated in Fig. [4] where cr q is calculated from Eq. ( I I 21 for 
HS cations and anions of varying size in TIP3P water 
in the absence of counterions 15311 . The main observa¬ 
tion here is that cr q significantly exceeds in the absolute 
magnitude the prediction of the continuum electrostat¬ 
ics. This means KQ«)I ^ I<71 in Eq. CS, which should 
lead to an overestimate of the number of adsorbed ions 
when the standard equations for the screening of free 
charge carriers in electrolytes are applied to the analysis 
of the mobility data & 

The overpolarization of the water dipoles attached 
to the surface ions might have significant implications 
for the interpretation of the mobility data. Figure [4] in- 
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dicates that the microscopic orientational order of the 
water dipoles next to a positive ion will significantly 
enhance its effective electrokinetic charge determined 
from the mobility measurement. Correspondingly, a 
negative adsorbed ion will appear more negative in the 
particle mobility. Therefore, the actual concentration of 
adsorbed ions can be significantly lower than estimated 
from mobility. This observation might help to explain 
the disagreement between the electrokinetic measure¬ 
ments p, 3JH and surface-sensitive spectroscopies 
HEUinnSi" regarding the concentration of the sur¬ 
face adsorbed ions. The actual extent of overpolar¬ 
ization requires more extensive simulations in realistic 
electrolytes. One also should not underestimate the po¬ 
tential effect of corrugation of any real water-oil inter¬ 
face 0, which will affect the average contact RDF G 
in Eq. (IT2l >. 

Figure [4] indicates that surface charge densities of 
large positive and negative ions with \q\ = 1 charge at 
the center are close in magnitude. However, this out¬ 
come might not hold for small ions adsorbed at the 
surface of a large particle. The product p\G is gener¬ 
ally asymmetric between cations and anions because of 
the asymmetry in the charge distribution of the water 
molecule [60]. In that case, the orientational order and 
the corresponding surface charge density will not com¬ 
pensate between the oppositely charged adsorbed ions, 
and a non-vanishing cr q will follow even at the total zero 
charge. The observable consequence of this asymme¬ 
try would be a shift between the iso-electric point of 
electrokinetic mobility and the point of zero charge, as 
reported for some systems 0. Overall, the main re¬ 
sult of the general formalism summarized by Eq. CD 
and simulations performed here is that the charge of free 
carriers and the effective electrokinetic charge incorpo¬ 
rating the interfacial dipolar order can be significantly 
different. 


4. Experimental testing and conclusions 

In conclusion, we have derived a simple equation 
[Eqs. ( ITU and CD] relating the effective charge of a hy¬ 
drated nanoparticle to the orientational order in the in¬ 
terface and the water density in contact with the solute. 
Both parameters carry asymmetry between positive and 
negative charges. Therefore, the surface charge density 
cr q induced by the positive and negative free carriers will 
not compensate and produce an overall nonzero value 
even when the total charge is zero. The electrokinetic 
charge can be substantially enhanced by the dipolar or¬ 
der in the interface and the theory predicts a non-zero 



a (A 3 ) 

Figure 5: The order parameter p\ (a) and the surface charge density 
(To (b) vs the polarizability of a polarizable HS solute, carrying the 
isotropic polarizability a and dipole moment mo = 5 D, in TIP3P wa¬ 
ter 1331 (error bars show the uncertainties of calculations). The hori¬ 
zontal arrow indicates the photoinduced alteration of the polarizability 
that inverts the mobility of the nanoparticle. 

effective charge when the interface is spontaneously po¬ 
larized in the absence of free charge carriers (charge- 
free mobility). The values of the surface charge density 
derived from simulations of uncharged nanometer-size 
solutes are consistent or exceed those typically reported 
from mobility measurements. 

Our development traces in spirit the well-established 
mechanism of electrophoretic mobility due to ions of 
electrolyte. Both ions and the dipoles of the solution 
surrounding the colloidal particle are neutral overall. 
However, there is an excess of ions within the shear sur¬ 
face of the particle, which results in the electrokinetic 
charge. Similarly, due to specifics of the divergent in¬ 
terfacial polarization, there is an inbalance in the bound 
charge between the polarized liquid inside and outside 
of the shear surface. The excess bound charge needs to 
be added to the excess free charge to establish the effec¬ 
tive electrokinetic charge. 

The derivation is performed here for a spherical so¬ 
lute, where the geometry of the interface produces a di¬ 
vergent radial polarization field. The model is not di¬ 
rectly extendable to fiat interfaces studied by simula¬ 
tions in the past QlUh. While the polarization field 
is clearly inhomogeneous next to a planar surface, it of¬ 
ten demonstrates positive and negative spikes 115^ ,39!], 
which can potentially compensate each other when the 
field is applied parallel to the interface to produce elec¬ 
trophoretic flow. The force ( F x ) along the plane of the 
surface (x-axis) writes 

(F x ) = E X S f p b (z)dz, (15) 
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where S is the surface area. If the density of the bound 
charge pb(z) integrates to zero, there is no net force. In 
this regard, the roughness of the interface, as suggested 
by Knecht et al 111511 . can provide the required condi¬ 
tions for a divergent polarization field which cannot be 
reduced to a one-dimensional integral shown above. 

The direct connection between the mobility of 
nanoparticles in water and the orientational order of the 
water dipoles in the interface offers opportunities for 
testing this prediction by experiment. One possible di¬ 
rection is the modification of the surface with chemi¬ 
cal groups (surface dipoles) altering the interfacial order 
[29]. We, however, recently discovered another prop¬ 
erty dramatically affecting the interfacial dipoles: the 
polarizability of the solute. Increasing the solute polar¬ 
izability drives the solute-water system to the point of 
instability of harmonic fluctuations expressed in terms 
of the solvent electric field inside the solute as the or¬ 
der parameter. Reaching the point of global instability 
toward fluctuations drives a structural transition of the 
hydration layer, which reorients the water di pole s and 
creates a high density of dangling OH bonds ]32[|. The 
emergent new structure of the interface also suggests, 
according to Eq. (fTO . the alteration of the sign of <t 0 


(q = 0). 

The results of MC simulations of HS solutes with 
changing isotropic dipolar polarizability a at the so¬ 
lute’s center are presented in Fig. [3 The size of the so¬ 
lute is maintained constant and only the polarizability is 
varied. One observes a switch from a positive to a neg¬ 
ative surface charge with increasing polarizability. In 
other words, the isoelectric point of electrophoretic mo¬ 
bility can be reached, and crossed, by manipulating the 
polarizability of the dissolved particle. This observation 
opens the door to experimental testing of the model. Po¬ 
larizability of semiconductor nanoparticles can be dra¬ 
matically increased by photoexcitation [33], which is 
predicted to invert the nanoparticles’ mobility (horizon¬ 
tal arrow in Fig. 
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Appendix A. Onsager reciprocal relations 

In order to prove the Onsager reciprocal relations 
for the problem of electro-osmotic current, one needs 
to consider the volume transport V in response to the 
applied gradient of the external electrostatic potential 
V</; cxt and the streaming current I in response to the ap¬ 
plied pressure gradient V/;: V = /-i 2 W) cxt , I = L 21 V/?. 
The Onsager reciprocal relations then require L \2 - L 21 ■ 


We start with the equation of motion for the stationary 
flow of an incompressible fluid (V • v = 0) along the z- 
axis of a capillary ^62] 

- 77 V 2 v z + piy • V)v z = -V z p + f z . (A.l) 


Here, v z (x,y) changes only along the cross section of 
the capillary ( x,y axes) and, therefore, V 2 = d 2 jdx 1 + 
d 2 /dy 2 . Further, 77 is the viscosity and p is the fluid den¬ 
sity. Since no convective motion of the liquid occurs, 
(v ■ V)v z vanishes. 

In contrast to the standard textbook description con¬ 
sidering free charges only, the force density f z in Eq. 
(ED is caused by a constant external field, E z = 
-V-(/> cxt , applied to the entire, free and bound, charge: 
f z = p{r)E z , p = pi - V • P. Since the curl of P dis¬ 
appears in the divergence V • P, one can put P = -Vc/>/, 
with the results 

P = --^-V 2 0, (/> = <pi- Amp b , (A.2) 

47T 

where 0 , is the electrostatic potential of free charges. 

We now proceed to calculate v z under the action of the 
force f z assuming no pressure applied to the capillary. 
The result from Eqs. ( IA. 1 b and (IA.2I) is 


V* 


~ 0 E 
Am) 


(A.3) 


where (j >0 is the potential at the shear surface at which 
v z = 0. In standard notations <po - e s £, where 'C is the 
^-potential and e s accounts for the screening by bound 
charges. Here, the potential of bound charges is a part 
of <f> and e s does not appear explicitly. A similar line of 
arguments can be applied to the potential of free charges 
<!>, connected to cf> through a closure relation. When the 
constitutive relations of continuous dielectrics are used, 
one has <p = e v 0 ,, where <!>, can be determined from solv¬ 
ing the Poisson-Boltzmann equation for the electrolyte 
next to the capillary wall. These details are irrelevant 
to our purpose since the derivation requires only the 
Coulomb law and the corresponding Laplace equation. 

From Eq. (IA.3I >. one gets the volume transport 


V = 



— Ti2V-0 ex t 


(A.4) 


with 

Ln = -^f(l-M<Po)dS. (A.5) 

We now turn to the streaming current when the cap¬ 
illary is subjected to the pressure gradient -V z p. The 
current is given by the equation 




-<f>)V 2 v z dS. 


(A. 6 ) 
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We now put f z = 0 in Eq. (IA. Ill , which results in 


I = L21V./7. 

(A. 7 ) 

It is easy to see that 


L 21 - L\ 2 , 

(A. 8 ) 

where L12 is given by Eq. ( 1 A . 5 b. 
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